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Abstract 

f-*) . The polarization tomography problem consists of recovering a matrix function 

OO ! / from the fundamental matrix of the equation Drj/dt = TT^frj known for every 

geodesic 7 of a given Riemannian metric. Here vr^ is the orthogonal projection 
■ onto the hyperplane 7-*-. The problem arises in optical tomography of slightly 

tf^ . anisotropic media. The local uniqueness theorem is proved: a C^-small function 

/ can be recovered from the data uniquely up to a natural obstruction. A partial 
global result is obtained in the case of the Euclidean metric on M'^. 
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■ 1 Introduction 



X 



First of all we shortly recall the physical motivation of the problem. See Section 5.1 of [5] 
for a detailed discussion. 

We consider propagation of time-harmonic electromagnetic waves of frequency a; in a 
medium with the zero conductivity, unit magnetic permeability, and the dielectric per- 
meability tensor of the form 

Eij = r?bi^ + -Xij, (1.1) 

where /c = cu/c is the wave number, c being the light velocity. Here n > is a function 
of a point x G M'^, and the tensor Xij = Xijip^^ determines a small anisotropy of the 
medium. The smallness is emphasized by the factor 1/fc. Equation (1.1) was suggested 
by Yu. Kravtsov [2J. By some physical arguments [3], the tensor x must be Hermitian, 

Xij — Xji- 
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In the scope of the zero approximation of geometric optics, propagation of electromag- 
netic waves in such media is described as follows. Exactly as in the background isotropic 
medium, light rays are geodesies of the Riemannian metric 

dt'^ ^n^{x)\dx\^; (1.2) 

the electric vector E(x) and magnetic vector H{x) are orthogonal to each other as well 
as to the ray; and the amplitude = \E\'^ = \H\'^ satisfies A = C/VnJ along a ray, 
where J is the geometric divergence and the constant C depends on the ray. The only 
difference between a slightly anisotropic medium and the background isotropic one consists 
of the wave polarization. The polarization vector rj — n~^A~^E satisfies the equation 
(generafized Rytov's law) 

Dri i ^. 

along a geodesic ray 7(t). Here t is the arc length of 7 in metric (1.2), 7 = dj/dt is the 
speed vector of 7, nj is the orthogonal projection onto the plane 7-*-, and D/dt — 7'^Vfe is 
the covariant derivative along 7 in metric (1.2). The right-hand side of (1.3) is understood 
as follows: vr^ and x ^-re considered as linear operators on T^^^-j = C^, and n^xV is the 
result of action of the operator ha^/x the complex vector G 7"*". Here 7-*- is the two- 
dimensional complex vector space consisting of complex vectors orthogonal to the real 
vector 7 e = T^(t) C T^^^^ = C^. Introducing the notation 

/ = ^X. (1.4) 

we rewrite (1.3) in the form 

^ = ^,fv. (1.5) 

Observe that / is a skew-Hermitian operator, fij = —fji- 

Let us now consider the inverse problem. Assume a medium under investigation to 
be contained in a bounded domain D (Z M.^ with a smooth boundary. The background 
isotropic medium is assumed to be known, i.e., metric (1.2) is given. The domain D 
is assumed to be convex with respect to the metric, i.e., for any two boundary points 
Xq, Xi G dD, there exists a unique geodesic 7 : [0, 1] — D such that 7(0) = xq, 7(1) = Xi. 
We consider the inverse problem of determining the anisotropic part Xij of the dielectric 
permeability tensor or, equivalently, of determining the tensor / on (1.5). To this end we 
can fulfill tomographic measurements of the following type. For any unit speed geodesic 
7 : [0,1] —>■ D between boundary points, we can choose an initial value 770 = ?7(0) G 7"'"(0) 
of the polarization vector and measure the final value r/i = ri{l) G 7^(/) of the solution to 
equation (1.5). In other words, we assume the linear operator 7^(0) — >• 7^(0) ^ Vi to 
be known for every unit speed geodesic 7 : [0, /] — > between boundary points. Instead 
of (1.5), we will consider the corresponding operator equation 

^ - Wit). (1.6) 

where f^^it) : 7"'"(t) 7"'"(^) is the restriction of the operator 7r^(t)/(7(t)) to the plane 
7-'-(i), and the solution is considered as a linear operator U{t) : 7"'"(t) — > 7"'"(0- Equation 
(1.6) has a unique solution satisfying the initial condition 

U{0) = E, (1.7) 
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where E is the identity operator. Since /^(j) is a skew- Hermit ian operator, the solution 
U{t) is a unitary operator. The final value of the solution 

mil) = m^GL{^^{i)) 

is the data for the inverse problem. Given the function $[/] on the set of unit speed 
geodesies between boundary points, we have to determine the tensor field / = {fij{x)) on 
the domain D. 

We consider the inverse problem in a more general setting. Instead of a domain D C 
with metric (1.2), we will consider a compact Riemannian manifold (M, g) of an arbitrary 
dimension n > 3, and an arbitrary complex tensor field / = (fij) on M. In such a setting, 
equation (1.6) makes sense along a geodesic 7. We will subordinate the manifold {M,g) 
to some conditions that guarantee smoothness of the data $[/] in the case of a smooth /. 

The two-dimensional case of = 2 is not interesting since and U become scalar 
functions and the solution to the scalar equation (1.6) is given by an explicit formula in 
this case. Therefore the inverse problem is reduced to the inversion of the ray transform 
I on second rank tensor fields, see the remark before Theorem 5.2.1 of [5J. 

Equation (1.6) can be slightly simplified. For a point x G M, let T^M be the com- 
plexification of the tangent space T^M. Instead of considering the operator U{t) on 7^(t), 
we define the linear operator 

Uit) : T^^(,)M ^ r^^(,)M 

by 

f/(t)i,x(,) = f/(t), u{mt) = i{t). 

If U{t) satisfies (1.6)-(1.7), then U{t) solves the initial value problem 

^ = (n^fn^)U, UiO)=E. (1.8) 

Equation (1.8) is more handy than (1.6) since all operators participating in (1.8) are 
defined on the whole of T^^^^M. The inverse problem consists of recovering the tensor 
field / from the data $[/](7) = [/(/) known for every unit speed geodesic 7 : [0, /] — > M 
between boundary points. Let us emphasize that the inverse problem is strongly nonlinear, 
i.e., the data $[/] depends on / in a nonlinear manner. 

The three-dimensional case, n = dimM = 3, is of the most importance for applications 
as we have shown above. On the other hand, the three-dimensional case is mathematically 
the exceptional one because, for a skew-symmetric /, the solution to the inverse problem 
is not unique. The non-uniqueness is discussed in Section 4. 

The main result of the present article is the local uniqueness theorem: the solution to 
the inverse problem is unique (up to a natural obstruction in the three-dimensional case) 
if the tensor field / is C^-small. See Theorem 5.1 below for the precise statement. 

Our method of investigating the inverse problem is a combination of approaches used 
in [9j and in Chapter 5 of [5j . First of all, following [9] , we reduce our nonlinear problem to 
a linear one as follows. Let /j (z = 1, 2) be two tensor fields and Ui{t) be the corresponding 
solutions to (1.8) with / = /j. Then u = U{^U2 — E satisfies 

Du 

-^=Pvr^(/2-/i)vr^g, u(0) = 0, (1.9) 
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where p = f/f^ and q = U-2. We consider p and q as operator-valued weights which are 
close to the unit operator if fi are C^-small. Assuming the weights p and q to be fixed, 
the solution u{t) to the initial value problem (1.9) depends linearly on / = /2 — /i. We 
study the linear inverse problem of recovering the tensor field / = /2 — /i from the data 
-^[/](7) = ""(0 given for all geodesies 7 : [0,/] — M between boundary points. In the 
case of a symmetric tensor field / and of unit weights, this linear problem was considered 
in Chapter 5 of [5]. We will demonstrate that the same approach works in the case of an 
arbitrary / and of weights close to the unit one. 

There is one more opportunity to extract a linear inverse problem from equation (1.8). 
Indeed, if W{t) = det U(t) is the Wronskian, then the function (f{t) = \nW{t) satisfies 

= tr (tt^/tt^). 

Therefore, for every unit speed geodesic 7 : [0, Z] — > M between boundary points, the 
integral 

I 

S[f]h) = J tr(vr^(,)/(7(t))7r^(,))rft (1.10) 


is expressed through the data $[/] by the formula 

^[/](7) = lndet$[/](7). 

The data S[f] depends linearly on /. Of course, some information is lost while the data 
$[/] is replaced with S[f]. In particular, S[f] is independent of the skew-symmetric part 
of/. 

We finally note that main results of the article are new and nontrivial in the case of 
M C M", n > 3, with the standard Euclidean metric. If a reader is not familiar with the 
tensor analysis machinery on the tangent bundle of a Riemannian manifold, he/she can 
first read the article for the latter simplest case. 

The article is organized as follows. Section 2 contains some preliminaries concerning 
Riemannian geometry and tensor analysis. In particular, we define some class of Rieman- 
nian manifolds for which the problem can be posed in the most natural way. Instead of 
considering the ordinary differential equation (1.8) along individual geodesies, we intro- 
duce a partial differential equation on the unit tangent bundle and pose an equivalent 
version of the problem in terms of the latter equation. In Section 3, we consider the 
corresponding linear problem and prove the uniqueness for weights sufficiently close to 
the unit in the case of n > 4. Section 4 discusses the three-dimensional case. In Section 
5, we check that the weights p and q are sufficiently close to the unit for a C^-small / and 
prove our main result. Theorem 5.1, on the local uniqueness in the nonlinear problem. 
In the final Section 6, we investigate the question: to which extent is a symmetric tensor 
field / determined by data (1.10). We give a complete answer to the question in the case 
of M = with the Euclidean metric. 
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2 Posing the problem and introducing 
some notations 



A smooth compact Riemannian manifold (M, g) with boundary is said to be a convex 
non-trapping manifold (CNTM briefly) if it satisfies two conditions: (1) the boundary 
dM is strictly convex, i.e., the second fundamental form 

II(e,0 = (V€^,0 for ^eT,{dM) 

is positive definite for every boundary point x G dM, where u is the outward unit normal 
vector to the boundary and is the covariant derivative in the direction ^; and (2) for 
every x G M and 7^ ^ G T^M, the maximal geodesic 'yx,i{t) determined by the initial 
conditions 7^,^(0) = x and ■jx,^{0) = ^ is defined on a finite segment [r_(x, ^), r+(x, ^)]. 
In what follows, we use the notations 7^;^^ and t±{x,^) many times. They are always 
understood in the sense of this definition. 

Remark. In [5] , the term CDRM (compact dissipative Riemannian manifold) is used 
instead of CNTM. In the case of M C M" with the standard Euclidean metric, this 
definition means that M is strictly convex. 

By TM = {(x, ^) I X G M, ^ G T^M} we denote the tangent bundle and by 

= G TM I lep = = 9^A^)ee = n, 

the unit sphere bundle. Its boundary can be represented as the union dflM = d^QM U 
d-ilM, where 

d±QM = {(x, e I a; G DM, ±(^, z/(x)) > 0} 

is the manifold of outward (inward) unit vectors. If 7 : [0, /] ^ M is a unit speed geodesic 
between boundary points, then (7(0), 7(0)) G 9_f2M and (7(/),7(Z)) G d+ilM. 

By T^M we denote the complexification of the tangent space T^M. The metric g 
determines the Hermitian scalar product on T^M 

{VX)=9^JVV■ (2.1) 

For ^ ^ G T,M, by Tj;^M = {r] e T^M \ (r/, ^) = 0} we denote the orthogonal 
complement of C, and by tt^ : T^M — ^ T^M, the orthogonal projection onto T^^M. 

Let tJM be the bundle of complex tensors that are r times contravariant an s times 
covariant. Elements of the section space C°°(rJM) are smooth tensor fields of rank 
(r, s) on M. In the domain of a local coordinate system, such a field u G C°°(rJM) 
can be represented by the family of smooth functions, u = the coordinates 

of u, where each index takes values from 1 to n = dimM. The metric g determines 
canonical isomorphisms tJM = Tq^^M = t^_^__^M. We will consider the isomorphisms as 
identifications. So, we do not distinct contra- and covariant tensors but use contra- and 
covariant coordinates of the same tensor. For example, for u G C°°{t^M) = C°^{tIM) = 

k- -k kl 

Uij = giku.j = gjkUi. = gikgjiu . 
In particular, such a tensor field determines the linear operator 

u{x) : T^M T^M, {ur,)' = u'^r]^ 
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at any point x G M. The product of two such operators is written in coordinates as 
{uv)ij = UikV^-. The dual operator has the coordinates u*^ = Uji. The operator is 
Hermitian (symmetric) if and only if Uij = Uji {uij = uji). The scalar product (2.1) 
is extended to tensors by the formula {u,v) = and determines the norm 

\u\^ = {u,u). For u,v E C°°{tIM), the norm of the product satisfies \uv\ < y/n\u\\v\, 
where n = dim M. 

We will also widely use semibasic tensor fields introduced in [3], see either Section 
3.4 of [5] or Section 2.5 of [7] for a detailed presentation. Let /3JM be the bundle of 
complex semibasic tensor fields of rank (r, s). It is a subbundle of Tg{TM) isomorphic to 
the induced bundle 7r*(rJM), where vr : TM — M is the projection of the tangent bundle. 
A tensor u G T^^,^^'^{TM) at (x, ^) G TM is semibasic if it is "pure contravariant in the 
^-variable and pure covariant in x" , i.e., 

■ ■ d d ■ 

For U C TM, by C°°{l3lM; U) we denote the space of smooth sections over U. The 
notation C°°(/5gM; TM) is abbreviated to C°°(/3^M). In the domain of a local coordinate 
system, such a field u G C°°{PIM) can be represented by the family of its coordinates, u = 
(iij^'[]j^{x, ^)), which are smooth functions of 2n variables (x, ^) = (x^, . . . , x", . . . , ^"). 
All the content of the previous paragraph is extended to semibasic tensor fields, where 
g remains the metric on M in the identification of contra- and covariant tensors. In 
particular, u G C°°{I3\M) determines the linear operator u{x,^) : T^M T^M for every 
(x,^) G TM. There are two important first order differential operators 

v,v: c~(/?:m) ^ c°°(/5:+,M) 

which are called the vertical and horizontal covariant derivatives. The operators are 
defined in local coordinates by the formulas 



d 

k "'3i-]s-> 



h . . d d 

r s 

I yia h-.-ia-ipia + l-.-ir _ "pp „.h--ir 

Z— ^ kp jl-js / J kja jl-ja-ipja + l-js^ 

a=l a=l 

where F*^ are Christoffel symbols. See Sections 3.4-3.6 of [5] for properties of these 

V h 

operators. Note that V= d/d^ and V= d/dx in the case of M C M" with the standard 
Euclidean metric and of Cartesian coordinates. 
The operator 

H = CVi : C^iPlM) C^iPlM) 

is of the most importance in the present article. It is called the differentiation with respect 
to the geodesic flow. 
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Given a tensor field / e C°°{t^M) on a CNTM {M,g), let us consider the boundary 
value problem 

HUix,0^nfi^)^^Uix,0 on QM, U\e_nM ^ E, (2.2) 

where E is the identity operator. A solution U = U{x, ^) is assumed to be a section of the 
bundle 0lM over QM, i.e., U e C{(3lM; QM). In the case of M C R" with the standard 
Euclidean metric, / and U can be considered as n x n-matrix valued functions oi x & M 
and of {x,$^) G M x §"^^ respectively. Problem (2.2) has a unique solution. Indeed, if we 
restrict (2.2) to an orbit of the geodesic flow, i.e., if we set x = 7(t) and ^ = 7(t) for a unit 
speed geodesic 7 : [0, /] — > M with 7(0) e dM, then we immediately arrive to the initial 
value problem (1.8). The boundary value problem (2.2) is thus equivalent to the family 
of initial value problems (1.8) considered for all unit speed geodesies simultaneously. The 
inverse problem is now formulated as follows: one has to recover the tensor field / given 
the trace 

m = U\9^nM (2.3) 

of the sohition to (2.2). 

In order to abbreviate further formulas, let us introduce the operator on tensors 
which maps f{x) to 7r^/(x)7r^ for (x,^) G QM, and write (2.2) in the shorter form 

HU = {P^f)U, U\9_nM = E. (2.4) 

Because of the factor and of the boundary condition on 9_QM, the solution U to (2.4) 
satisfies 

U{x,0^-U%x,m^C (2.5) 
Therefore the non-trivial part of the data (2.3) consists of the restrictions $[/](a;, ^)|yx m 
for (x,^) G d^QM. This agrees with the above discussion of the relationship between 
(1.6) and (1.8). The solution U is continuous on QM and C°°-smooth on QM \ Q[dM) 
as one can easily prove using the strict convexity of the boundary. 

Concluding the section, let us mention one more inverse problem that is not considered 
in the present article. Let gl{T^M) be the space of all linear operators on T^M. The 
operator P^ participating in (2.4) is the orthogonal projection of the space gl{T^M) onto 
the subspace 

{/ G gl{T^M) I /e = n = 0}. 
Let us introduce the smaller subspace 

{/ G gl{T^M) I /e = n = 0, tr/ = /; = 0} 

and denote by the orthogonal projection onto the latter subspace. The corresponding 
inverse problem for the equation 

HU={Q^f)U, U\o_nM = E (2.6) 

is also of a great applied interest. To explain the physical meaning of (2.6), let us return 
to equation (1.5) considered in the three-dimensional case for a skew-Hermitian tensor 
/. The polarization vector rj on (1.5) is a complex two-dimensional vector subordinate to 
one real condition I77I = 1. Therefore 77 can be described by three real parameters. Two of 
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these parameters can be chosen to determine the shape and position of the polarization 
elhpse on the plane 7"*", while the last parameter is the phase of the electromagnetic 
wave. See section 6.1 of [5] for a detailed discussion of the subject. Only the first two 
of these parameters are measured in practice. Deleting the wave phase from the data is 
mathematically equivalent to replacing the operator with Q^. The authors intend to 
consider the corresponding inverse problem for (2.6) in a subsequent paper. 

3 Linear problem 

Let (M, g) be a CNTM. Choose two semibasic tensor fields p,qe C°°{(3lM; VtM) satisfy- 
ing 

P*(a;,0^ = e, g(^,0^ = e (3.1) 

For a tensor field / G C°°(r^^M), consider the boundary value problem on VtM 

Hu=p{P^f)q, u\9_nM = 0. (3.2) 

The problem has a unique solution u G C{PlM; QM) and, in virtue of (3.1), the solution 
satisfies 

u{x,O^ = u*{x,O^ = 0. (3.3) 

In this section, we consider the inverse problem of recovering the tensor field / from the 

data 

F[f] = ulo^nM. (3.4) 

The factors p and q on (3.2) are considered as weights. We will assume the weights to 
be close to the unit weight E in the following sense: the inequalities 

V V 

\p-E\<e, \q-E\<e, \Vp\ < £, \Vq\ < £ (3.5) 

hold uniformly on QM with the norm | ■ | defined in Section 2. The value of e will be 
specified later. 

Equation (3.2) is initially considered on QM. To get some freedom in treating the 
equation, we extend it to the manifold T^M = {{x,^) G TM | ^ 7^ 0} of nonzero vectors. 
The weights are assumed to be positively homogeneous of zero degree in ^ 

p{x,t^) = p{x,^), q{x,t^) = q{x,^) for t > 0. 

Then the right-hand side of (3.2) is positively homogeneous in ^ of zero degree because / 
is independent of ^. The solution u must be extended to T°M as a homogeneous function 
of degree —1 

u{x, to = t'^U{x, for t>0 

because the operator H increase the degree of homogeneity by 1. 

Let us discuss smoothness properties of the solution u. It can be expresses by the 
explicit formula 
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p(7(t), 7(t)) Aw/(7(t))g(7(t), 7(i)) dt, 



where 7 = 71,.,^ and T^'° is the parallel transport of tensors along the geodesic 7 from 
the point 7(t) to 7(0) = x. The integrand is a smooth function. Therefore smoothness 
properties of u are determined by that of the integration limit r_(x, ^). The latter function 
is C°°-smooth on T°M \ T{dM) but has singularities on T^{dM). Therefore some of 
integrals considered below are improper and we have to verify their convergence. The 
verification is performed in the same way as in Section 4.6 of [S]. So, in order to simplify 
the presentation, we will pay no attention to these singularities. 

Besides (3.5), we will impose some smallness condition on the curvature of (M, g). For 
(x, ^) e fiM, let K{x,S,) be the supremum of the absolute values of sectional curvatures 
at the point x over all two-dimensional subspaces of T^M containing ^. Define 



k{M,g)= sup j tK{j,^^{t),%^^{t))dt. (3.6) 



Theorem 3.1 For any n > A, there exist positive numbers 6 = 6{n) and e = e{n) such 
that, for any n- dimensional CNTM (M, g) satisfying 

k{M,g)<6 (3.7) 

and for any weights p, q G C°°{PlM; flM) satisfying (3.1) and (3.5), every tensor field f G 
C'^{tIM) can he uniquely recovered from the trace (3.4) of the solution to the boundary 
value problem (3.2) and the stability estimate 

\\f\W<C\\F[f]U (3.8) 

holds with a constant C independent of f . In the case of n = 3, the same statement is 
true for a symmetric tensor field f . 

In the case of a real symmetric / and unit weights, this theorem is a partial case of 
Theorem 5.2.2 of [5j. We will show that the same proof works with some modifications 
for Theorem 3.1. 

Proof of Theorem 3.1. We rewrite equation (3.2) in the form 

Hu = P^f + r, (3.9) 

where 

r = {p- E)P^f + pP^f{q - E). (3.10) 

The remainder r is small by (3.5). Because of (3.3), the function u = {uij{x,C,)) is 
orthogonal to ^ in both indices 

CUiJ=eU^j = 0. (3.11) 

We write down the Pestov identity for the semibasic tensor field u (see Lemma 4.4.1 
of [3] for the case of a real u and Lemma 5.1 of [6] for the general case) 

2Re{\/u,VHu) = \Vu\^ + \/iv' + Viw' -TZilu], (3.12) 

where (■, ■) and | ■ | are the scalar product and norm on semibasic tensors defined in 
Section 2, 

v' = Re i^eV'u'''' ■ VjUi.i, - f VV^'^ ■ Vi^ni,) , (3.13) 
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= Re W''' • ^j^hi^ , (3-14) 

n,[u] = Rkpi.e^'^V'^u^''' ■ + Re {{R;',juP'' + RlljU^^n^^W'ui,,,) , (3.15) 

and (Rijki) is the curvature tensor. Identity (3.12) holds on any open D C TM for every 
u e C\PlM;D). In our case, D = T^M \ T{dM) and e C~(/3jM;D) as is shown 
above. 

The most part of the proof deals with the left-hand side of (3.12). We will first 
transform it by distinguishing some divergent terms and then will estimate it. 
Prom (3.9) 

(W, VHu) = {Vu, V{Pif)) + {Vu, Vr). (3.16) 

We will first investigate the first term on the right-hand side of (3.16). To this end we 
represent / as 

where (fij) is a semibasic tensor field orthogonal to ^ in both indices 

hC = he = 0, (3.18) 
semibasic covector fields a and b are orthogonal to ^ 

aie^biC = 0, (3.19) 

and c{x,C,) is a scalar function. One can easily check the existence and uniqueness of the 
representation. The (vector versions of the) fields a and b are expressed through / by the 
formulas 

b=j^^7rJX (3.20) 

As follows from (3.17)-(3.19), 

Pd = f 

or in coordinates 

= fij = fij - «iO - ^j^ - Ciiij. (3.21) 

Differentiating the last equality with respect to ^ and using the fact that / is independent 
of ^, we obtain 

V V V - V - 

Vk{Pe,f)ij = -^jVu^i - e'^kbj - e^jVkC - gjkcii - Qikbj - {gik^j + gjk^A 
Therefore 

{Vu,V{Pif)) = VV^' ■ VkiPjhj = 

h , . . / V V V \ 

= V u^^i- ^jVkai - e'^kbj - e^jVkC - Qjkai - Qikbj - {gik^j + gjke)c\ ■ 

h .. 

The tensor V is orthogonal to ^ in the indices i and j as follows from (3.11). Therefore 
the last formula is simplified to the following one: 

{Vu, V{Pd)) = -V^Uip ■ a' - V%i • b\ 
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h h 

Introducing the semibasic covector fields Siu and by the equahties 

i5iu)i = V^Uip, i52u)i = V^'upi, (3.22) 
we write the result in the form 



{Vu, V(n/)) = -{6iu, a) - {62U, b). (3.23) 



This implies the estimate 

2 Re{Vu, V(Ps/)) < ^{\6iu\' + \62u\') + + \b\'), (3.24) 

where (3 is an arbitrary positive number. 

h 

Next, we transform the expression \Siu\ by distinguishing a divergent term 



tq 



= Vpiu'^V'ui,) - u'PVpV'u,, = Vpiu'^V'u,,) - ulVpVqu''. (3.25) 
By the commutator formula for horizontal derivatives (see Theorem 3.5.2 of [5]), 



Substituting this value into the previous formula, we obtain 

iSiMp = Re ( - ufVqVpu''' + Vpiu'^^V^UiS) + 7^2M, (3.26) 



where 



7^2[n] = Re (ul{R%^eVku"' - R)pqU'' - R 
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We now transform the first summand on the right-hand side of (3.26) in the order 
reverse to that used in (3.25) 

\5iu\^ = Re ( - Vgiul^P^'") + ■ VpW"^ + Vpiu'^Wig)) + 7^2M. 

Introducing the semibasic vector field vi by the formula 

Re (u^'V%k - UjkV''u^'^ , (3.27) 



we write the result in the form 



iSiMp = Re (vV^ ■ VkUj.^ + ViiviY + 7^2M. (3.28) 



In the same way, we obtain 

h „ / h . ., h \ h 



\52u\^ = Re (y'u^" ■ Vju.k) + Viiv2y + TZsiu] (3.29) 
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with 



and 



{V2y = Re (u'^v'^Ukj - Ukj^v'^u'^^ (3.30) 



7^3M = Re (<(it:J,^e v.f^^^ - R%,u^' - Ri,,u'^')) . 

Taking the sum of (3.28) and (3.29), we have 

\Siu\'^ + \S2u\'^ = Rc (^VV'' • VjUik + VV'' • Vit%) + Viv' + n4[u\, (3.31) 

where 

V = Vi+ V2, 7^4[-u] = 7^2N + T^3[u]. 
Introduce the semibasic tensor field z = (zijk) by the formula 

^iUjk = -^iHu)jk + Zijk- (3.32) 

The idea of this new notation is that the tensor z is orthogonal to ^ in all its indices 

^ Zijk — Zijk — ^ Zijk — 0, (3.33) 

h h 

while the tensor 'S/u — (ViUjk) has the mentioned property only in the last two indices. 
The summands on the right-hand side of (3.32) are orthogonal to each other, so 

\Vu\^ ^ j^\Hu\^+\z\^. (3.34) 

The first two terms on the right-hand side of (3.31) can be expressed through z. Indeed, 
one easily see with the help of (3.32) and (3.33) that 

Re (y'u^'' ■ VjUik + V'u^^ ■ VkUj^ = Re {z'^^Zjik + z'^^z^ji) <2\z\'^. 
With the help of the last inequality, (3.31) implies the estimate 

iSiMp + \52u\'^ < 2\zf + \J,v' + 7^4M (3.35) 



which, together with (3.24), gives 

2 Re(V«, V(P^/)) < /^kr + |(|ar + \h\^) + ^4^^ + ^-Tl^- 

Substitute values (3.20) for a and h into the last formula 

2Re(W, V(P^/)) < M + + Wn\^) + f + f ^4N. (3.36) 

Next, we estimate the second term on the right-hand side of (3.16). We differentiate 
equality (3.10) with respect to ^ taking the independence / of ^ into account 

Vr = V ((p - E)Pi) f + V{pP^)f{q - E) +pPJVq. 

12 



In what follows in the proof, we denote different constants depending only on n = dim M 
by the same letter C. On using (3.5), we obtain from the last formula 

Ce 

|Vr|<^|/|. 



From this 



2Re(VM, Vr) < Ce{\Vu\^ + ^l/l')- 



Combining (3.36) and (3.37), we obtain from (3.16) 



(3.37) 



(3.38) 

Estimating the left-hand side of the Pestov identity (3.12) by (3.38), we obtain for 
lel = l 

|Vwp+V.^^-/3|zp-|(|7r5/ep+|7r^reP)-C^^(|Vup + |/P) < V.(§^^-^^0+7^N, (3.39) 



where 



n[u] = 7^l[M] + ^TZiiu]. 



We multiply inequality (3.39) by the volume form dT, = \dxUj{C,) A dV"{x)\, integrate over 
flM, and transform the integrals of divergent terms by Gauss-Ostrogradskii formulas (see 
Theorem 3.6.3 of [5]) 



nM 



+ _ 2)\Hu\' -/3\z\'- |(kc/er + k^rer) - Cei\Vu\' + \f\') 
< j {^v-v,u)dT?'^-^ + J n[u]dE. 



(3.40) 



dQM QM 

The second term on the left-hand side has appeared because w is positively homogeneous 
of degree —1 in ^ and satisfies {^,w) = as is seen from (3.14). Substituting the 

value |V^P = + from (3.34), we write the result in the form 



< J {^v-v,u)dT?'''^ + j 7^[^t]dS. 



dT. < 



(3.41) 



ami fiM 

Assuming P < I, integrals on the right-hand side of (3.41) can be estimated exactly 
as in Section 5.5 of [5] 



n\u]dE 



7lM 



<Ck{M,g) / \Vu\^dT<C6 / IV^I^ciS 



(3.42) 
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< D\\u\ 



a+nM\\m, 



(3.43) 



where k{M,g) is defined by (3.6). The second inequality on (3.42) is valid because of 
(3.7). The constant D on (3.43) depends on {M,g), unlike the constant C on (3.42) 
which depends only on n. 



Combining (3.41)-(3.43) and using again the equality |Vm| 



we obtain 



I ]^(^i_p_Ce-CS)\z\' + {n-l-Ce-C5)\Huf-^iMC\' + \n^rC\')-C^^^^ 

QM 



— -D||M|g^f7M||/i-i 



(3.44) 



with some new constant C depending only on n. 

Let us compare \Hu\ and \P^f\- The estimate |r| < Ce\f \ follows from (3.5) and 
(3.10). The latter, together with (3.9), imphes 



\Hu\'>\P^f\'-Ce\f\'. 
Using the last inequality, we transform (3.44) to the final form 

{l-f3-Ce- C5)\z\^ + {n-l-Ce- C5)\PJ[' 



(3.45) 



Let us remind that /3 is an arbitrary number satisfying < /3 < 1. 



(3.46) 



Lemma 3.2 For every Riemannian manifold (M, g) of dimension n> A and every point 
X e M, the Hermitian form 



B 



(/,/)= / 



i)ip^/i^-2(ivei^+i7r^rer)H^.(o 



is positive definite on the space of second rank tensors at x. Moreover, the estimate 

B(f,f)>c\f\' 

holds with a positive constant c depending only on n. In the case of n = 3, the same 
statement is valid for symmetric tensors. 

The proof of the lemma will be given later, and now we finish the proof of Theorem 
3.1 by making use of the lemma. 
By the lemma, the inequality 



cii ri|i2 = c j \f\^dv^(x) < 

M M n^M 



[1 - P - Ce - CS)\z\'^+ 
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+ ^n-l-Ce- C5)\P^f\^ - + WFi?) - Ce\ff\ dcu^iOdV^ix) (3.47) 

holds for P = 1, e = 6 = 0. By continuity and by estimates \P(f \ < |/|, Itt^/^I < 
I/I; k?/*'^! < I/I foi' 1^1 = 1) the same inequality holds for some positive c,e,6, and /3 
independent of / and satisfying 1 — P — Ce — C5 > 0. Combining (3.46) and (3.47), we 
obtain 

\\f\\h < C\\u\a,nM\\l. = C\\F\f]\\l. 
with C = D /c. This finishes the proof of Theorem 3.1. 

Proof of Lemma 3.2. One easily checks the equality -B(/, /) = B{u, u) + B{v, v) for 
a complex tensor / represented as / = u + if with real u and v. Therefore it suffices to 
prove the statement for a real tensor /. 

The corresponding symmetric bilinear form is 

B{f,h)= j [{n-l){P^f,h)-2{{Tx^fiM) + {^iri.h*i))\du^{i). 

Obviously, B{f, h) = for a symmetric / and skew-symmetric h. Therefore, it suffices to 
prove the positiveness of B on the spaces of real symmetric and skew-symmetric tensors 
separately. 

The positiveness of B on the space of real symmetric tensors is proved in Lemma 5.6.1 
of [5], where vr^/^ is denoted by P^j^f- 

Thus, we have to consider the quadratic form B on the space of real skew-symmetric 
tensors. On making use of an orthonormal basis, we identify T^M with M" endowed with 
the standard scalar product and identify Q^M with the unit sphere fl C M". So 

-S(/,/) = i / \in-l)\PJ\'-A\nJ^\']dco{0, (3.48) 
UJ u) J I J 

n 

where uj is the volume of f2 and duj is the standard volume form on f2. 
For a skew-symmetric / 

Hfi = /e (3.49) 

since is orthogonal to ^. Therefore formulas (3.20) and (3.21) are simplified to the 
following ones: 

a = -b = /e, {P^f)^, = - (/04- + Ufi)i 

for 1^1 = 1. Using the last equality and {fC,,0 = 0; we easily calculate 

\PJ\' = \f\'-m\' for |el = l. (3.50) 
Substitute (3.49) and (3.50) into (3.48) 

ii?(/,/) = (n-l)|/r-^^^ / m'dco. (3.51) 

UJ J 

On using the obvious relation 

/ C.C . rl,., — } 

1/n for i = j, 
15 



00 

h 



we find 



i/i/«r^.=ii/r. 

Inserting this value into (3.51), we see that 



iB(/./) = !?!^^|/|l 
CO n 



This implies the positiveness of B on skew-symmetric tensors for n > 4. 



4 Three-dimensional case 

We will first show that both our problems, linear and nonlinear, possess some non- 
uniqueness in the three-dimensional case. 

Let (M, g) be a three-dimensional CNTM which is assumed to be oriented. Every tan- 
gent space T^M is a three-dimensional oriented Euchdean space. So, the vector product 

T^M X T^M T^M, {v,w)^vxw 

is well defined. It is extended to the C-bilinear operation 

T^M X T^M T^M, {v,w)^vxw 

on complex vectors. For a complex vector field v e C°°{tqM), we denote by Ly e 
C°°{t^M) the operator of vector multiphcation by v, 

Ly{x)r] = v{x) XT] for rj e T^M. 

Note that is a skew-symmetric tensor field. Quite similarly, for a semibasic vector field 
V e C°°{PqM), the operator e C°°{/3lM) is defined. Let us prove the formula 

P^Ly = n^L^n^ = ^L^. (4.1) 

We remind that (•, ■) and | • | are defined in Section 2. This is a pure algebraic local 
formula. So, we can use a positive orthonormal basis (ci, 62, 63 = C/ICI) i^i T^M. In such 
a basis 

1 \ f -V3 V2 \ / -1^1 \ 





/ 1 














1 





















-V3 


V2 









-V2 


Vl 







and the formula follows immediately. 
Next, we prove the formula 

VL^ = 0. (4.2) 

h 

The formula is quite expectable since V{ = 0. Nevertheless, it needs a proof. Here, we 

h 

have to use a general coordinate system since V is a differential operator. The vector 
product is expressed by the formula 

[V X wy = -^{Vi+iWi+2 - Vi+2Wi+i) 

v9 
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in general coordinates, where g — det(gijj) and indices are reduced modulo 3. Therefore 

(L,),,,+i = -^v'+' (4.3) 

and (4.2) is proved by straightforward calculations with making use of ViQjk — and 
Vi^' = 0. 

We derive from (4.1)-(4.2) and the formula H{l/\^\'^) = that 



for a complex function A G C°^(M). Thus 

H{^h) = P^Ly,. (4.4) 

This gives us the following non-uniqueness in the linear problem. If the function A vanishes 
on the boundary, X\dM — 0, then u{x,^) — AL^/|^p solves the boundary value problem 

Hu = P^f, Ud_nM = (4.5) 

with / = LvA and satisfies 

F[f] = ula^nM = 0. (4.6) 

The boundary value problem (4.5) coincides with (3.2) for the unit weights p — q — E, 
and (4.6) means that / cannot be recovered from F[f]. 

The same arguments give us some non-uniqueness in the nonlinear problem. Fix 
a function A e C°°{M) vanishing on the boundary, X\qm — 0, and, for (x, ^) e JIM, 
define the linear operator U{x,^) on T^M whose matrix in a positive orthonormal basis 
(61,62,63 = ^) is 

cos X{x) — sin X{x) 
sinA(x) cosA(a:;) 
1 

For a real function A, U{x, ^) is the rotation of T^M around the axis ^ by the angle A(a;). 
In the case of a complex A, the operator is also well defined although its geometric sense is 
more comphcated. The semibasic tensor field U e C°°{I3\M] QM) satisfies the equation 

HU = {P^L^x)U (4.7) 

and boundary condition 

U\dnM = E. (4.8) 

Indeed, (4.8) is obvious. Let us prove (4.7). In virtue of (4.1), equation (4.7) is equivalent 
to the following one: 

HU{x,0 = {^Kx).i)LiU{x,i). (4.9) 

Let 7 be a unit speed geodesic. Setting x = j{t), ^ — ^{t) in (4.9), we arrive to the 
equation 

DU(t) dX(t) 



dt dt 



Ht)U{t), (4.10) 
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where A(t) = A(7(t)) and U(t) = U{'j{t),'y(t)). Conversely, if (4.10) holds for any unit 
speed geodesic 7, then (4.9) is true. To prove (4.10), we choose an orthonormal basis 
{ei(t) , e2(t) , e^it) = ■j(t)) of T^(^t)M which is parallel along 7. In such a basis, (4.10) is 
equivalent to the matrix equation 



cosA(t) -sinA(t) 0\ d\ f ^ ~^ ^ \ / cosA(t) -sinA(t) 
sinA(t) cosA(t) =— 1 sinA(t) cosA(t) 
l/ \000/VO 1 



which is obviously true. 

Comparing (4.7)-(4.8) with (2.3)-(2.4), we see that a field / = Lva with X\qm = 
cannot be recovered from data (2.3). 

Let us introduce the definition: / G C°°{tIM) is said to be a potential field if it can 
be represented as / = Lva for some function A G C°°{M) vanishing on the boundary, 
MdM = 0. Potential fields constitute the natural obstruction for the uniqueness in the 
both, linear and nonlinear, problems. We are going to prove that a solution to the linear 
problem is unique up to the obstruction. The corresponding local result for the nonlinear 
problem will be obtained in the next section. First of all we prove 

Lemma 4.1 (on decomposition). Let {M,g) be a compact oriented three-dimensional 
Riemannian manifold with boundary. Every tensor field f G C°°{tIM) can be uniquely 
represented as 

f = Lyx + f (4.11) 

with some A G C°°(M) satisfying 

\\dM = (4.12) 

and some tensor field f G C°°{tIM) satisfying the condition: the 2-form fijdx^ A dx^ is 
closed, 

difijdx' A dx^) = 0. (4.13) 

The summands of decomposition (4.11) are called the potential and closed parts of / 
respectively. Note that (4.13) involves only the skew-symmetric part of /, i.e., a symmetric 
tensor field is closed. Lemma 4.1 can be derived from the Hodge- Morrey decomposition 
[8] but we give a shorter independent proof. 

Proof of Lemma 4.1. We first prove the uniqueness statement. Assume (4.11)- 
(4.13) to be valid. Applying the exterior derivative d to the form fijdx^ A dx^ and using 
(4.11) and (4.13), we obtain 

d{^Lyx)ijdx^ A dx^^ = d{fijdx^ A dx^). 

On using (4.3), one can check by a straightforward calculation in coordinates that 

d(^{Ls/x)ijdx' A dx^^ = -2A\^/gdx^ A dx'^ A dx^, 

where A is the Laplace-Beltrami operator. Thus, the function A solves the Dirichlet 
problem 

AA = -^(§% + §% + ff), AU = 0, (4.14) 

V OX^ OX^ OX'^ / 
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where /~ is the skew-symmetric part of /, i.e., — ^{fij — fji)- The solution to the 
Dirichlct problem is imique. 

The existence statement is proved by reverse arguments. Given /, let A be the solution 
to the Dirichlet problem (4.14) and / = / - Lva- Then (4.11)-(4.13) holds. The lemma 
is proved. 

We restrict ourselves to considering the inverse problem for closed fields only. 

Theorem 4.2 There exist such positive numbers 6 and e that, for any oriented 3-dimen- 
sional CNTM {M, g) satisfying (3.7) and for any weights p,q e C°°{(3lM;flM) satisfying 
(3.1) and (3.5), every closed tensor field f e C°°{tIM) can he uniquely recovered from the 
trace (3.4) of the solution to the boundary value problem (3.2) and the stability estimate 

WfWh < C mf]\\l. + WfUh^ ■ \\F[f]h^) (4.15) 
holds with a constant C independent of f . 

Proof follows the same line as the proof of Theorem 3.1 with the following difference: 
the left-hand side of the Pestov identity (3.12) will be estimated in a different way by 
making use of the closeness of /. 

We represent / as the sum of symmetric and skew-symmetric fields 

f = f+ + f- f.+ = f + f .T = -f- 

Taking the symmetry of Christoffel symbols into account, the closeness condition for / 
can be written as 

V./-. + V,/,l + V./^ = 0. (4.16) 

The vector f~^ is orthogonal to ^ and therefore 7T^f~^ = f~^. Formulas (3.20) take now 
the form 

« = + ^/"^' ^ = ^^^^"""^ - j^^r^- (4-17) 

We write equation (3.2) in the form (3.9) with the remainder r defined by (3.10). Then 
we write the Pestov identity (3.12) for u with terms defined by (3.13)-(3.15). The left 
hand-side of the identity can be written as in (3.16). 

h V 

The main problem is estimating the first term (Vu, V(-P^/)) on the right-hand side of 
(3.16). To this end we represent / in form (3.17) with a and b defined by (4.17). Then 
(3.23) holds. In view of (4.17), equation (3.23) can be written as 

{Wu,V{PJ)) = -j^{Siu + 52U,7rJ+C) - j^{Siu-52U,f-0- (4-18) 
We transform the second term on the right-hand side of (4.18) by distinguishing di- 

h h 

vergent terms. Using definition (3.22) of Siu and S2U, we write 

{Siu - S2U, f-0 = iy.u'" - v.unfik^' = 
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Now, we use the closeness condition (4.16) to transform this formula to the following one: 
Finally, we distinguish a divergent term from the first summand on the right-hand side 



This can be written as 

h h _ h 



{5iu - 52U, ro = {Hu, n + m'^'^fki + c'u^'n + eu'^jz). (4.19) 



Next, wc calculate the first term on the right-hand side of (4.19) by making use of 
(3.9) and (3.21) 

{Hu, n = {HuY'iij = mff + ni-, = (r - a^e - - cce + ni^- (4.20) 

Since /~ is a skew- symmetric tensor, c^^^-^ f~j = 0. This means that the term c^*^-' on 
the right-hand side of (4.20) can be omitted. By the same reason, the term can be 
replaced with (/")*■', the vector a = j^n^f^ can be replaced with j^fC, and the vector 
b can be replaced with j^f*^. In such the way, (4.20) takes the form 

{Hu, n = ir r - ^(/%e + /'^e^e)/^- + /-)• 

The second term on the right-hand side is independent of f'^, and the formula can be 
written as 

{Hu,r)^\r?-^,\re + {r,n- 

Substitute (4.21) into (4.19) 

{5,u - '52U, ro = i/-r - + (r, /-) + ^(evV"^, + eu^^ij, + euv-,,). 

Inserting this expression into (4.18) and estimating the first term on the right-hand side 
of (4.18) in a similar way as in deriving (3.24), we arrive to the inequality 

2Re(VH,V(P,/)) < ^{{Lu\' + \S2u\') + ^j7:,re- 

' + ^irer + 4(^3)^ + ^|r|- I/I, (4.22) 

2 



\cr ' ' ' 1^1 

where 

(^^3)^ = MCu^'fk, + eu^'f]k + i'^'^fW (4-23) 

and (3 is an arbitrary positive number. The last term on the right-hand side of (4.22) can 
be estimated by C£|/p/|^p as follows from (3.5) and (3.10). Estimating the first term 
on the right-hand side of (4.22) by (3.35), we obtain 

2Re(W,V(n/))</3|^r+^|7r,/+er-^l/1'+^l/-eP+|^l^ 

(4.24) 
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with the same curvature dependent term TZ^lu] as in (3.36) and 



V = 2^'"! +'"2) +"^3, 

where Vi {i = 1,2,3) are defined by (3.27), (3.30), and (4.23) respectively. 

The second term on the right-hand side of (3.16) is estimated by (3.37) as before. 
Combining (3.37) and (4.24), we obtain from (3.16) 

2Re{Vu,VHu)<P\z\' + ^M+^f-2\r\' + Mr^\' + Ce{\ff+\^^^^^ 

(4.25) 

for 1^1 = 1. 

We estimate the left-hand side of the Pestov identity (3.12) by (4.25) and write the 
result in the form 

(l-Ce)\Vu\'+ViW^-P\z\'~\n^re+'^\r^ 

Integrating this inequality and transforming the integrals of divergent terms in the same 
way as in (3.40), we obtain 



[1 - Ce)\Vu\' + \Hu\' -P\z\'- || + 2|/-r - 41/-^!' - Ce\f\' 



< 



Substituting the value |VMp = + from (3.34), we write the result in the form 

J [(1 - /3 - Ce)\z\' + (2 - Ce)\Hu\' - + 2\r\' - 41/-^^ - Ce\ff]dJ: < 



(4.26) 



The curvature dependent integral on (4.26) is estimated as before in (3.42) 



7^l[^^]c^E 



<C5 j \Vu\'di:, 

QM 



(4.27) 



while the boundary integral is estimated in a little bit different way. Namely, instead of 
(3.43), we have the estimate 



{v — V, i/)dE 



2n-2 



mM 



< D(^\\u\oj^fiM\\m 



(4.28) 
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The second term on the right-hand side of (4.28) appears because of the dependence of 
■5 on / as is seen from (4.23). Inequahty (4.28) is proved in the same way as estimate 
(4.7.2) of 15]. 

Combining (4.26)-(4.28) and using again the equahty | Vwp = \z\^ + \Hu\'^, we obtain 



(l_^_C'5-C'e)|;.p + (2-C5-C'5)|iftzp-l|V+eP+2|rp-4|rep-Ce|/p' 



Q.M 



< D[\\u\Q^nM\\]ii + ||/|aAf||L2 • Ikb+f^MlUa)- (4.29) 

The tensors and /~ are orthogonal to each other as well as P^^f^ and P^f^ are 
orthogonal to each other. Therefore 

\f? = \n' + \f-\', in/P = inn' + inr 

With the help of these equalities, (3.45) gives 

\Hu\^ > - Ce\f+\^ - Ce\f-\\ (4.30) 

We use (4.30) to transform the estimate (4.29) to the final form 

^l_P_C5-Ce) J \z\'dJ:+ J \^(2-C5-CE)\P^r\'-^M^^\'-Ce\f- 



+2 J (|/1'-2|rep-C'£|rnrfs<D(||M|a+nA/||^i + ||/|aA/||L-||«|a+nA/||L2) (4.31) 

QM 

with some new constant C. 

For /3 = 1 and e = 6 = 0, the left-hand side of (4.31) is the integral over M of the 
Hermitian form 

A{fj) = B{r,n+2Q{r,f-), 

where B is the same as in Lemma 3.2 and 

Q{fJ)= I i\f\'-m\')dco,{0- (4.32) 

By Lemma 3.2, the form B is positively definite on symmetric tensors in the 3-dimensional 
case. The form Q is positively definite on skew-symmetric tensors. Indeed, repeating the 
arguments of the proof of Lemma 3.2, we derive in the 3-dimensional case 

- / {\f\'-m\')du;^{0 = \\f\' 
uj J 3 

for a skew-symmetric tensor /. Now, the proof is finished in the same way as in Theorem 
3.1. 
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5 Nonlinear problem 

We return to considering the inverse problem of recovering a tensor field / G C°°{tIM) 
on a CNTM (M, g) from the data 

m = Ula^nu, (5.1) 
where U G C{j3lM; QM) is the solution to the boundary value problem on QM 

HU = {PJ)U, U\a_nM = E. (5.2) 
We will prove the uniqueness under the following smallness assumptions on /: 

\f{x)\<e for xg9M, (5.3) 



j |/(7,.,5(t))|rft<£, j |V/(7.,«(t))Mt <£ for (x,0g5+^^M. (5.4) 

T-{x,i) T-{X,S,) 

We remind that we use notations introduced in Section 2. In particular, V is the covariant 
derivative. Note that these smallness conditions are quite similar to that of Theorem 2 of 

m- 

Theorem 5.1 It is possible to choose a positive number 6 = S{n) for n > 4 such that, 
for an n- dimensional CNTM {M,g) satisfying the curvature condition (3.7), there exists 
a positive number e = e{M,g) such that the following statement is true. Let two tensor 
fields fi G C°°{tIM) {i = 1,2) satisfy (5.3)-(5.4) and $j = <l>[/j] be the corresponding 
data. Then the estimate 

||/2-/i||l2 <C7||$-i$2-^||hi (5.5) 

holds with a constant C independent of fi. In particular, fi = f2 if = $2- In the case 
of n = 3, the same statement is true under the additional assumption that f2 — fi is a 
closed tensor field. 

Proof. Let Ui G C{(3lM:,^lM) [i = 1,2) be the solution to the boundary value 
problem 

HU, = iPJ,)U,, mo_nM = E. (5.6) 
According to (5.2), the solution satisfies 

= = e (5.7) 

Set u = f/f - E. Using the equalities HE = 0, H{U{^U2) = HU^^ ■ U2 + U^^HU2, 
and HU^^ = -U^^HUi ■ U{\ one easily derive from (5.6) that u solves the boundary 
value problem 

Hu = U{\P^if2-fi))U2, u\9_nM = (5.8) 

and 

u\d+nM = $r^'^2 - E. (5.9) 



23 



Setting / = /2 — /i, P = ^1 ^, 1 = U2, we write (5.8)-(5.9) in the form 

Hu=piPJ)q, u\9_nM = 0, (5.10) 

F[f] = ula+nM = $r''^2 - E. (5.11) 

We have arrived to the hnear problem considered in Sections 3 and 4. If we will prove 
that the weights p = f/f ^ and q = U2 satisfy conditions (3.1) and (3.5), we would be able 
to apply Theorems 3.1 and 4.2 to obtain the statement of Theorem 5.1. Condition (3.1) 
is satisfied by (5.7). The weight p = solves the boundary value problem 

Hp = -piP^fi), p\d-nAi = E 

which is very similar to (5.2). Theorem 5.1 is thus reduced to the following 

Lemma 5.2 Let {M,g) be a CNTM and a tensor field f G C^^tIM) satisfy (5.3)-(54). 
Then the solution U to the boundary value problem (5.2) satisfies the estimates 

\U-E\<Ce, \VU\<C£ on VlM (5.12) 

with some constant C depending on {M,g) but not on f. 

To prove Lemma 5.2, we need the following estimate for solutions to linear ordinary 
differential equations (see Lemma 4.1 of Chapter IV of [Ij): 

Lemma 5.3 Let y = {yi{t), . . . ,yN{t)) be the solution to the initial value problem 

^ = Ait)y + fit), yiO) = yo, 

where f{t) is an N -dimensional vector and A{t) is an N x N -matrix. Then 

t t 

\y{t)\ < (bo| + j |/(r)|rfr)exp j |A(r)|rfr, 


where \ A{t)\ is the operator norm of the matrix A{t) defined with the help of the standard 
norm \- \ on C^. 

Let us adjust this statement to our geometric setting. 

Lemma 5.4 Let {M,g) be a CNTM and f G C°°(/5^M; fiM), A G C^{Pl^M;nM) be 
two semibasic tensor fields. Let y G C(/3^M;f2M) be the solution to the boundary value 
problem 

Hy = Ay + f, y\d_nM = yo 
with some yo G C{l3^M; d-VtM). Then 



Il/(x,0l<c(|yo(7(r-(^,0),7(r-(a^,0))l+ / |/(7(t), 7^)1^^) x 

T-{X,S,) 



xexp[c j \A{^{t),^{t))\dt 

for (x,^) G QM, where 7 = •jx,^- Here the norm \ ■ \ is defined in Section 2, and the 
constant C depends on m and {M,g). 
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The constant C appears in Lemma 5.4 since different norms arc Tiscd in this Lemma 
and Lemma 5.3. In what follows in this section, we denote different constants depending 
on (M, g) by the same letter C. 

Proof of Lemma 5.2. Let U be the solution to (5.2). Then U — E solves the 
boundary value problem 

H{U -E) = {PJ){U -E) + PJ, {U - E)\9_aM = 0. 

Applying Lemma 5.4, we obtain the estimate 



\{U-E){x,0\<C J \P^^t)f{l{t))\dt e^p[c J \P^^t)f{l{mdt\< 

T- {X,£,) T- {X,0 





<C j Ifi^mdtexp [C J \fm)\dt\, 



T- {X,£,) T- {X,(,) 

Together with (5.4), this implies the first of inequalities (5.12). 

V h 

The proof of the second of estimates (5.12) is more troublesome because 'SJU and VU 

h 

must be estimated together but Vf^ is unbounded near Vl{dM). 

h 

We start with estimating V^^ on d-VtM. To this end we consider equation (5.2) at 
a boundary point (x, ^) e 9_QM. Because of the boundary condition U\q_^m = E, 
equation (5.2) gives 

{HU)\9_aM^Pif. (5.13) 

Let us choose boundary normal coordinates (a;^,...,x") in a neighborhood V of the 
boundary point such that the boundary is determined by the equation a;" = 0, > in 
V , and Qin = Sin- Because of the boundary condition Uj\xn=o = Sj, the equalities 

Va^j|a:"=O = (0 < « < 71 - 1) 

hold, and equation (5.13) becomes 

Since = —{^, where u — ^{x) is the unit outward normal to the boundary, this gives 
with the help of (5.3) 

m(x,0\<j^^ for (x, e 9_QM, (i., 7^ 0. (5.14) 

V 

yU vanishes on d^QM as follows from the condition U\q_^m — E, 

VU\9_nM = 0. (5.15) 

V h 

Applying the operators V and V to equation (5.2), we obtain 

WHU = {P^f)VU + {WP^fU, VHU = (Pe/)VC/ + (P«V/)C/. (5.16) 
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V h 

We have used that V/ = and V/ = V/ since / is independent of ^. 
Operators V and H satisfy the commutator formula 

VH = H\/ + \/ (5.17) 

. h V h 

as follows immediately from the definition H = ,^*Vi and from the fact that V and V 

h 

commute. The commutator formula for V and if is a little bit more complicated. Indeed, 
using the commutator formula for horizontal derivatives (Theorem 3.5.2 of [5j), we see 
that 

h ■ h h ■ h h 



This can be written as 

VHU = HWU + TZi[WU]+n2[U] (5.18) 

with some algebraic operators TZi and 7^2 on semibasic tensors which are determined by 
the curvature tensor. The operator TZ2 satisfies TZ2[E\ = 0. Therefore the first of estimates 
(5.12), which is already proved, implies the inequality 

\n2[U]\ = \n2[U -E]\<Ce (5.19) 

with some constant C depending on the curvature bound. 

Using commutator formulas (5.17) and (5.18), we write (5.16) as 

H{VU) = {PJ)VU -VU + F, (5.20) 

H{\/U) = -7^l[Vf/] + (Pc/)W + G, (5.21) 

where 

F={wP^)fU, G={P{s/f)U -n2[U]. (5.22) 
We first consider (5.20)-(5.21) as a linear system of ordinary differential equations in 

V h 

coordinates of '^U and '^U with free terms F and G. Applying Lemma 5.4 to the system, 
we obtain the estimate 

|Vf/(a:,OI < C||vf/(7(r-(x,0),7(r-(x,0))l + |W(7(r-(a;,0),7(r-(x,0))l + 



+ j (|F(7(t),7W)l+|G(7(t),7W)l)rft}exp [c j {\n^{{^{t),^m+\Pmfiim+\E\)dt 

T-{x,i) r-(x,^) 

(5.23) 

The first summand of the expression in braces is equal to zero by (5.15). In virtue of 
(5.14), the second summand of this expression is estimated as 

|Vf/(7(r-(^,0),7(r-(^,0))l< 



|(7(r_(x,0),K7(^-(a;,0))r 
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By Lemma 4.1.2 of [S], the estimate 



|(7(r_(x,0),K7(r-(^,0)))l 
holds. Combining two last estimates, we obtain 



< C 



|Vf/(7(r-(^,0),7(r-(^,0))l< 



r-{x,0\' 

As is seen from (5.22), (5.19), and (5.4), the integral inside the braces in (5.23) can be 
estimated by Ce. Finally, the integral under the exponent on (5.23) is estimated by some 
constant. Therefore (5.23) implies the estimate 

m{x,0\ < I for {x,OenM. (5.24) 

Now, we consider (5.20) as a linear system of ordinary differential equations in coordi- 

V h 

nates of with the free term — + F. Applying Lemma 5.4 to the system and using 
the homogeneous initial condition (5.15), we obtain the estimate 



|W(x,OI<C j (|W(7(t),7(t))l + |i^(7(t),7(t))l)t^texp [c J \P^^t)fiimdt 

T_(x,$) T_(a;,{) 

(5.25) 

In virtue of (5.24), the first integral on the right-hand side of (5.25) can be estimated by 
Ce. Estimating then the second integral by (5.4), we obtain the second of inequalities 
(5.12). The lemma is proved. 



6 Kernel of the operator S 

In this section, we restrict ourselves to considering symmetric matrix functions on the 
whole of endowed with the standard scalar product (■,■). 

Let M(3) be the space of complex-valued symmetric 3 x 3-matrices. Such a matrix 
/ e M(3) is considered as the linear operator / : — > C^. By §^ we denote the unit 
sphere in and by TS^ = {(x,0 | ^ e S^,x G M^, (x,^) = 0}, the tangent bundle of 
the sphere. Given ^ G let ^-^ = {i] E \ = 0} be the complex two-dimensional 

space of vectors orthogonal to C, and vr^ : ^ C'^, the orthogonal projection onto C,'^. 

Let iS(M^;M(3)) be the Schwartz space of M(3)-valued functions on R^. The linear 
operator 

S : 5(M^;M(3)) C°°(r§2) 

is defined by 

oo 

S[f]{x,0= j trMi^ + t0^i)dt for (x,0eT§2, (6.1) 

— oo 

compare with (1.10). We are going to answer the question: to which extent is a symmetric 
matrix function / G iS(M^;M(3)) determined by the data S[f]7 Since S[f] depends 
linearly on /, the question is equivalent to studying the kernel of the operator 5*. 
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Let us recall the definition of the ray transform 

oo oo 

nfK^^o = j {fix+m,odt= jux+toeedt for {x,oeTs^. (6.2) 

—00 —00 

See Chapter 2 of [5J for the theory of the ray transform on the Euclidean space. Our 
approach to studying the kernel of S is based on the following observation: S[f] = — /[/] 
for a sjTiimetric matrix function / with zero trace, see equation (6.7) below. 

Theorem 6.1 A symmetric matrix function f = {fij{x)) G 5(R^;M(3)) belongs to the 
kernel of the operator S if and only if it satisfies the system of partial differential equations 

Rllf] :=/ll;ll + 2/l2;12 + /22;22 + /sS;!! + /33;22 " | 

R2[f] :=/ll;ll + 2/l3;13 + /22;11 + /22;33 + /33;33 — (6-3) 
Rslf] '■= /ll;22 + /ll;33 + 2/23;23 + /22;22 + /33;33 = J 

where the tensor notations for partial derivatives fij-ki = d"^ fij/dxkdxi are used for brevity. 

The theorem gives the full system | i = 1,2,3} of local linear functionals that 

can be recovered from the data S[f]. 
Proof. Represent the matrix / as 

f = f + \E, tr/ = 0, (6.4) 

where E is the unit matrix and A G iS(]R^) is a scalar function. Then 

S[f] = S[f] + S[XE]. (6.5) 

As follows from definition (6.1) of the operator S, 

S[XE] = 2I[XE], (6.6) 

where 

oo 

I[XE]{x,0= j X{x + tC)dt for (a;,0er§2 

— oo 

is the ray transform of the matrix function XE. 
If (^, ?7, Q is an orthonormal basis of M'^, then 

= tr / = (/(x)e, + (/>)^, V) + (/(^)C, C) 

and 

tr(7r^/>)7r5) = (/>)//,//) + (/>)C,C). 

This implies 

tr(7r5/> + tOvr5) = ~{hx + ti%i). 
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Integrating the last equality with respect to t and recalling definitions (6.1) and (6.2) of 
S and /, we see that 

S[f] = -/[/]. (6.7) 
Substituting (6.6)-(6.7) into (6.5), we obtain 

S[f] = -I[f-2XE]. 

Thus, / is in the kernel of S if and only if 

I[f-2XE] = 0. (6.8) 

Now, we apply Theorem 2.2.1 of [5j which states that (6.8) is equivalent to the existence 
of a vector field v such that 

f-2\E = dv, (6.9) 



where 

d : C°°(M^M') -> C°°(M^M(3)), (dv) 



2 \dxj dxi 



is the inner derivative. Theorem 2.2.1 of j5] is formulated and proved for compactly 
supported tensor fields. Nevertheless, the same proof works for / — 2XE G 5(]R^; M(3)) 
and gives the vector field v belonging to the Schwartz space 5(M^; M.^). 

Let us express A and / through /. Applying the trace operator to the first of equations 
(6.4) and taking the second one into account, we see that 

A = ^tr/. (6.10) 

From (6.4) and (6.10), 

fij = fij - ^tif ■ 6ij. (6.11) 

Substitute (6.10)-(6.11) into (6.9) to obtain 

dv = f - ti f ■ E. (6.12) 

Equation (6.12) represents the overdetermined system of six first order partial differ- 
ential equations in three unknowns (fi, ^2, t's). The solvability condition for the system 
is presented by Theorem 2.2.2 of [S]: equation (6.12) is solvable if and only if the right- 
hand side of (6.12) belongs to the kernel of the Saint- Venant operator. Here, we prefer to 
use the version R of the Saint- Venant operator which is defined by the equation before 
formula (2.4.6) of [5]. So, the solvability condition for (6.12) is 

Rh = 0, where h = f -tif ■ E. (6.13) 

The operator R is defined by the formula 

A{Rh)ijki = hik-ji — hjk-ii — hii-jk + hji-ik- 

It possesses the following symmetries: 

{Rh)ijki = —{Rh)jiki = —{Rh)ijik = {Rh)kiij. 
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Because of the symmetries, the tensor Rh has six hnearly independent components, and 
equation (6.13) is equivalent to the system 



-Mf] 

-R4[f] 
-RrAf] 
-Rdf] 

Substitute the value hij = fij — (/n + /22 + fz^j^ij into the last system 



[Kn)i2i2 


— ^11;22 


— 2/ll2;12 


+ ^22;11 = 0, 




(i?/i)l313 


= ^11;33 


— 2/li3;i3 


+ hs-u = 0, 




(i?/l)2323 


= ^22;33 


— 2/;,23;23 


+ /i33;22 = 0, 




{Rh)l2l3 


= hn;23 


~ /il2;13 " 


- ^13;12 + ^23;11 


= 0, 


(i?/l)2123 


= ^22;13 


~ ^12;23 " 


~ ^23;12 + ^13;22 


-0, 


(i?/i)l323 


— ^33; 12 


— ^13;23 " 


~ ^23;13 + ^12;33 


= 0. 



Rlif] 

R2[f] 

Rsif] 

R4[f] 

Rdf] 

Reif] 



/ll;ll + 2/i2;12 + /22;22 + /33;11 + /33;22 = 0, 
/ll;ll + 2/l3;l3 + /22;11 + /22;33 + /33;33 — 0, 
/ll;22 + /ll;33 + 2/23;23 + /22;22 + /33;33 = 0, 
/l2;13 + /l3;12 + /22;23 ~ /23;11 + /33;23 — 0, 
/ll;13 + /l2;23 ~ /l3;22 + /23;12 + /33;13 — 0, 
/ll;12 — /l2;33 + /l3;23 + /22;12 + /23;13 = 0. 



(6.14) 



For a symmetric matrix function / = {fij{x)) G iS(]R^; M(3)), system (6.14) contains 
only three independent equations. More precisely: each of the last three equations of 
(6.14) can be obtained from the first three equations by taking linear combinations, dif- 
ferentiation, and integration. To prove this, let us rewrite system (6.14) in terms of the 
Fourier transform g{^) — f. Applying the Fourier transform to each equation of (6.14), 
we arrive to the system 



Ri[9] 

R2[g] 
Rsia] 

Re[g] 



ilgii + 2^16^712 + ei^?22 + {il + ei)^733 = 0, 



(6.15) 



ihii + 26e3^?i3 + {il + il)g22 + elfe = 0, 
(el + ^3)^11 + ilg22 + 266^23 + ihz^ = 0, 

iii^gu + Ci69'i3 + 665'22 - Cife + S^to = 0, 

iii^gii + 665'i2 - Cls'is + ^16523 + 665'33 = 0, 
:= iliigii - Cigi2 + 669'i3 + ei6fl'22 + 669'23 = 0. 

One can easily see the following three relations between equations of system (6.15): 

2664b] = ClRM + CiMg] - 6'4b], 

2664m = esRilg] - e2R2[g] + 6'4M, 
2664 b] = -esRiig] + e2R2[g] + 6'4b]- 

Therefore three last equations of (6.15) follow from three first equations at least if g{^) 
depends continuously on 6 
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Deleting three last equations from system (6.15), we obtain the equivalent system 



Ri[9] 

Ma] 



= ilgii + 2ei6^?i2 + el^?22 + (e? + el)^?33 = 0, 

= ihll + 266^13 + {il + ^3)^22 + ^1^33 = 0, !^ (6.16) 
= (CI + ^3)^11 + Cl^22 + 266^23 + ^3^33 = 0. 



The same is true for system (6.14): deleting three last equations from (6.14), we will 
obtain the equivalent system (6.3). The theorem is proved. 

System (6.16) enables us to answer the question: which integral moments of / can be 
determined from S[f\l Indeed, the Tailor series of the function g{E,) = f is 



^.•^(0 ~ E E b^^lifr. (6.17) 

where 



m=0 \a\=m 



\a\ = m 



are the integral moments of order m. Assuming / to be in the kernel of S, let us insert 
series (6.17) into system (6.16). Since the coefficients of (6.16) are homogeneous functions 
of ^, the system does not mix moments of different orders. This means that we can take 
g{$,) in the form 

■m 

\a\=m 

if we are looking for moments of order m. 

Let us start with considering zero order moments. We substitute the expressions 
Qij — jif^ into (6.16). Equating coefficients at the same degrees of ^ at the resulting 

equations, we easily find that nf^ — for every (i,j). This means that the integral 
/iR3 f{^)dx can be determined from the data S[f]. 

Next, we consider first order moments. We substitute the expressions 

9ij{o = + /^lj,k2 + /^lj,k3 

into system (6.16). Equating coefficients at the same degrees of ^ at the resulting equa- 
tions, we arrive to the system 

1^11,1 + 1^22,1 — 0) A*ll^,l + 1^33,1 — 0' 1^11,2 + 1^22,2 — 0) 

A*22^,2 + 1^33,2 — 0) A*ll'',3 + A'33^3 ~ 0' A*22'',3 + 1^33,3 ~ 0' 
1^11,2 + 2/i^2^,l + 1^33,2 — 0) A*n'',3 + '^1^13,1 + A*22^,3 ~ 0' f^il,3 + A*22^,3 + '^1^23,2 — 0) 

2^12^,2 + f^i"!,! + ^33^1 ~ 0, 2/i^3^^3 + fJ^22,l + f^33,l ~ 0' I^U,2 + '^1^23,3 + A''33^,2 ~ 0- 

The general solution to the system looks as follows: 

(1) _ (1) _ (1) _ (1) _ (1) _ n (1) _ 

«1) /^12,l — «2, /^13,1 — «3, /^22,1 — /^23,1 — U, /^33,1 — 
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Hll 2 — ~0'2i 1^12,2 ~ '^1' 1^1^,2 ~ 1^22,2 ~ '^2; ^^23, 2 ~ '^3; A''33,2 ~ ~'^2^ 

(1) _ (1) _ n (1) _ (1) _ (1) _ (1) _ 

/^ll,3 — ~'^35 /^12,3 — ^13 3 — fli, /i22,3 " ~'^3) A^23,3 — ^2, A^33,3 " '^3? 

where (01,02,03) are arbitrary constants. Eliminating the constants, we obtain the fol- 
lowing independent system of 15 linear combinations of first order moments of / which 
can be recovered from the data S[f]: 

where Sij is the Kronecker tensor. System (6.18) is considered for such {i,j,k) that at 
least two of these indices are different. A similar consideration is possible for integral 
moments jJ^^J^^^lf] of an arbitrary order m. 
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